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Abstract. We simulate cooling of superfluid neutron stars with nucleon cores where direct Urea process is forbid- 
den. We adopt density dependent critical temperatures T cp (p) and T cn (p) of singlet-state proton and triplet-state 
neutron pairing in a stellar core and consider a strong proton pairing (with maximum T™ ax £ 5 x 10 9 K) and 
a moderate neutron pairing (T™™ ~ 6 x 10 s K). When the internal stellar temperature T falls below T c ™ ax , the 
neutrino luminosity Lqp due to Cooper pairing of neutrons behaves oc T 8 , just as that produced by modified Urea 
process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing 
neutrino emission acts like an enhanced cooling agent. By tuning the density dependence T cn (p) we can explain 
observations of cooling isolated neutron stars in the scenario in which direct Urea process or similar process in 
kaon/pion condensed or quark matter are absent. 
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1. Introduction 

Thanks to Chandra and XMM-Newton missions, there is a 
great progress in observations of thermal radiation emer- 
gent from the surfaces of isolated (cooling) middle-aged 
neutron stars (e.g., Pavlov & Zavlin 2003). A comparison 
of these data with theoretical models of cooling neutron 
stars gives a method to constrain (still poorly known) fun- 
damental properties of supranuclear matter in neutron- 
star cores, such as the composition and equation of state 
of the matter and its superfluid properties. 

So far, the observations can be explained by a num- 
ber of vastly different theoretical models (e.g., Page 
UWhial I1998bl Tsuruta et al. VMM Khodelet al. IStTMl 
Blaschke et al. EDOl Yakovlev & Pethick l2TjOH and refer- 
ences therein). Particularly, one can employ the simplest 
models of neutron stars with the cores composed of nu- 
cleons (or nucleons/hyperons), or containing pion conden- 
sates, kaon condensates or quarks. The simplest model of 
a non-superfluid nucleon core which cools via modified 
Urea process of neutrino emission (without any power- 
ful direct Urea process) cannot explain the observations: 
some neutron stars (e.g., PSR B1055-52) turn out to be 
much warmer, while others (e.g., the Vela pulsar) are much 
colder than those given by this model. Warmer stars can 
be explained (Kaminker et al. 12001(1 assuming a strong 
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proton superfluidity in the core: such a superfluidity sup- 
presses modified Urea process and slows down the cool- 
ing. However colder stars require some cooling mechanism 
which is faster than the modified Urea process. 

Explanations of observations of colder stars presented 
in the literature invoke usually either a powerful direct 
Urea process in nucleon (or nucleon/hyperon) matter 
or similar processes in kaon-condensed, pion-condensed, 
quark matter in the inner cores of massive neutron stars. 

In this paper we present a new scenario of neutron star 
cooling. We adopt the simplest model equation of state 
of supranuclear matter in neutron star cores (Douchin 
& Haensel I20()l|l involving only nucleons, electrons and 
muons. This equation of state forbids direct Urea process 
in all stable neutron stars. We will show that the enhanced 
cooling required to explain colder isolated neutron stars 
can be produced by neutrino emission due to a moderately 
strong triplet-state pairing of neutrons. This new interpre- 
tation is possible only for a specific density dependence of 
the critical temperature of neutron pairing. 

In the next section we outline the observational basis; 
the cooling scenario is given afterwards. 

2. Observations 

Table 1 summarize observations of isolated (cooling) 
middle-aged (10 3 ^ t ^ 10 6 yr) neutron stars, whose ther- 
mal surface radiation has been detected (or constrained). 
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We present the estimated stellar ages t and effective sur- 
face temperatures T s °° (as detected by a distant observer). 

Two young objects, RX J0822-4300 and IE 1207.4- 
5209, are radio-quiet neutron stars in supernova remnants; 
RX J1856.4-3754 and RX J0720.4-3125 are also radio- 
quiet neutron stars. Other objects — the Crab and the 
Vela pulsars, PSR B1706-44, PSR J0538+2817, Geminga, 
and PSR B1055-52 — are observed as radio pulsars. 

RX J0205+6449 and the Crab pulsar are associated 
with historical supernovae and their ages are certain. For 
RX J0822-4300, we take the age of the host supernova 
remnant, Puppis A. As can be deduced, e.g., from a dis- 
cussion in Arendt et al. fl991), its age ranges from 2 to 5 
kyr; the central value is t = 3.7 kyr (Winkler et al. I1988|l . 
For IE 1207.4-5209, we also adopt the age of the host su- 
pernova remnant (G296.5+10). According to Roger et al. 
I|1988jl . it is t ~ 3 - 20 kyr. For the Vela pulsar, we take 
the age interval from the standard characteristic spindown 
age of the pulsar to the characteristic age corrected due 
to the pulsar glitching behaviour (Lyne et al. 1996). The 
age of PSR J0538+2817, t = (30 ± 4) kyr, was estimated 
by Kramer et al. (2003) from the measurements of the 
pulsar proper motion relative to the center of the host su- 
pernova remnant, S147. The age of RX J1856. 4-3754 has 
been revised recently by Walter & Lattimer l|2002|l from 
the kinematical reasons. Following these authors we take 
the central value t = 500 kyr and choose such an error-bar 
of t to clearly distinguish the revised value from the value 
t = 900 kyr reported previously by Walter H2()01(l on the 
basis of less accurate parallax measurement. The charac- 
teristic age of RX J0720. 4-3125 has been estimated by 
Zane et al. I|2002|l . Kaplan et al. H2002|l and Cropper et al. 
( 2004) from X-ray measurements of the neutron-star spin- 
down rate. We adopt the central value t = 1300 kyr with 
an uncertainty by a factor of 2. The ages of three other 
pulsars, PSR B1706-44, Geminga, and PSR B1055-52, are 
the characteristic pulsar ages assuming an uncertainty by 
a factor of 2. 

For two youngest sources, RX J0205+6449 and the 
Crab pulsar, no thermal emission has been detected, but 
the upper limits on surface temperature T s °° have been 
established (Slane et a!. IM)2l Weisskopf et al. l2lM|l . The 
surface temperatures of the next five sources, RX J0822- 
4300, IE 1207.4-5209, Vela, PSR B1706-44, and PSR 
J0538+2817, have been obtained using hydrogen atmo- 
sphere models (see references in Table 1). Such models are 
more consistent with other information on these sources 
(e.g., Pavlov et al. 2002J) than the blackbody model. On 
the contrary, for the Geminga and PSR B1055-52 we 
present the values of T s °° inferred using the blackbody 
spectrum because this spectrum is more consistent for 
these sources. 

Let us notice that from Table 1 we have excluded PSR 
B0656+14 which was considered earlier (e.g., Yakovlev et 
al. E002). A combined analysis of new X-ray and optical 
observations of the source (with the improved distance 
from new parallax measurements of Brisken et al. 2003) 
leads either to unrealistically small values of the neutron 



star radius (in the blackbody model) or to unreasonably 
small distance to the star (in the hydrogen atmosphere 
model); see, e.g., Zavlin & Pavlov (2002). This makes cur- 
rent interpretations of the data unreliable. 

The surface temperature of RX J1856. 4-3754 is still 
rather uncertain. A wide scatter of T s °°, obtained by dif- 
ferent authors, takes place because X-ray and optical ob- 
servations are not described by one blackbody model. This 
can be explained, for instance, by the presence of hot spots 
on the neutron star surface. Thus, we adopt the upper 
limit T s °° < 0.65 MK, which agrees with the value of T s °° 
obtained either with the "Si-ash" atmosphere model of 
Pons et al. (2002) or with the model of condensed sur- 
face layers of Burwitz et al. (2003). It agrees also with 
the model of nonuniform surface temperature distribution 
suggested by Pavlov & Zavlin ( 2003J). In the latter case, 
the mean surface temperature T s °° « 0.5 MK is below our 
upper limit of T s °°. 

Finally, T s °° for RX J0720.4-3125 is taken from Motch 
et al. H2U03[) who have interpreted the observed spectrum 
with a model of a hydrogen atmosphere of finite depth. 

For PSR J0538-4300, PSR B1055-52, and RX J0720.4- 
3125, the authors cited in Table 1 have not reported any 
error bars of T s °° . We adopt 20% uncertainties which seem 
to be appropriate for these sources. 

3. Physics input and calculations 

We will simulate cooling of neutron stars using our gen- 
erally relativistic cooling code described by Gnedin et al. 
( 2000) . We adopt a moderately stiff equation of state of 
neutron star interiors proposed by Douchin & Hacnsel 
(200l|. According to this equation of state, neutron star 
cores (regions of the densities p > 1.3 x 10 14 g cm~ 3 ) 
consist of neutrons, with the admixture of protons, elec- 
trons and muons. All constituents exist everywhere in a 
core, except for muons which appear at p > 2.03 x 10 14 
g cm -3 . The most massive stable star has the (gravita- 
tional) mass M = M max = 2.05 Af©, the central density 
p c = 2.9 x 10 15 g cm -3 , and the (circumferential) radius 
R = 9.99 km. The central densities and masses of eight 
neutron star models (with M from 1.111 Af© to 1.994 Af©) 
are presented in the right panel of Fig. 

All physics input is standard. The effects of muons are 
included as described by Bejger et al. ( 2003). We assume 
no envelope of light elements on stellar surfaces (Sect. EJ. 
The code calculates the cooling curves, which give the de- 
pendence of the effective surface stellar temperature T s °° 
on stellar age t. Let us remind that neutron stars are born 
hot in supernova explosions (with internal temperatures 
T ~ 10 11 K) but gradually cool down via neutrino emis- 
sion from the entire stellar body and via heat diffusion 
to the surface and thermal surface emission of photons. 
Qualitatively, one can distinguish three cooling stages. At 
the first ('non-isothermal') stage [t 100 yr) the main 
cooling mechanism is neutrino emission but the stellar in- 
terior stays highly non-isothermal. At the second ('neu- 
trino') stage (10 2 & t & 10 5 yr) the cooling goes mainly 
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Source 


t [kyr] 


T s [MKJ 


Confid. 


References 


PSR J0205+6449 


0.82 


<1.1 b) 


- 


Slane et al. ©02 1 


Crab 


1 


<2.0 b) 


99.7% 


Weisskopf et al. (I5U0H 


RX J0822-4300 


2-5 


1.6-1.9 a) 


90% 


Zavlin et al. ifTMty 


IE 1207.4-5209 


3-20 


1.4-1.9 a) 


90% 


Zavlin et al. 


Vela 


11-25 


0.65-0.71 a) 


68% 


Pavlov et al. (I5UDH 


PSR B1706-44 


~17 


n 09+O.OI a) 
U - OZ -0.34 


68% 


McGowan et al. (1200411 


PSR J0538+2817 


30 ±4 


~ 0.87 a) 




Zavlin & Pavlov lIMi^ 


Geminga 


~340 


~ 0.5 b) 


90% 


Zavlin & Pavlov lIMi^ 


RX J1856.4-3754 


~500 


<0.65 




see text 


PSR B1055-52 


~540 


-0.75 b) 




Pavlov & Zavlin 


RX J0720.4-3125 


~ 1300 


~ 0.51 o) 




Motch et al. (2003) 



Inferred using a hydrogen atmosphere model 
Inferred using the black-body spectrum 



via neutrino emission from isothermal interiors. At the 
third ('photon') stage (t 10 5 yr) a star cools predomi- 
nantly through the surface photon emission. 

The new element of our present studies is the equation 
of state of Douchin & Haenscl (2001). We have chosen it 
because it forbids the powerful direct Urea process of neu- 
trino emission (Lattimer et al. I1991[l in all stable neutron 
stars (M < M max ). In this case, a non-superfluid neutron 
star of any mass Mq < M < M max will have almost the 
same (universal) cooling curve T s °° (t) (the dotted curve in 
the right panel of Fig. QJ. At the neutrino cooling stage, 
this curve is determined by the neutrino emission due to 
the modified Urea process. The curve is almost indepen- 
dent of the equation of state of neutron star cores (Page 
& Applegate I1992[l as long as the direct Urea process is 
forbidden. As has been indicated by many authors (see, 
e.g., Yakovlev & Pethick EDM and references therein) and 
seen from Fig. ^ this universal cooling model is certainly 
unable to explain the data. For instance, it gives T s °° much 
lower than that of PSR B1055-52, but much higher than 
that of the Vela pulsar. We will show that all the data can 
be explained assuming superfluidity of neutron-star cores. 

It is well known that neutrons and protons in stel- 
lar cores can be in superfluid state. Proton superfluidity 
is caused by singlet-state proton pairing, while neutron 
superfluidity is produced by triplet-state neutron pair- 
ing. These superfluidities can be specified by density de- 
pendent critical temperatures for protons and neutrons, 
T cp (p) and T cn (p). Results of calculations of these tem- 
peratures from microscopic theories show a large scatter 
of critical temperatures depending on a nucleon-nucleon 
interaction model and a many-body theory employed. In 
particular, recently Schwenk & Friman l|2()()4|) and Zuo et 
al. (2004) have obtained weak neutron and proton pairing 
in neutron star cores but many other calculations give 
much stronger superfluidity (e.g., Lombardo & Schulze 



200l| also see references in Yakovlev et al. I1999fl . In this 
situation it is reasonable to consider T cp (p) and T cn (p) as 
unknown functions of p (consistent with predictions of mi- 
croscopic theories) which can hopefully be constrained by 
comparing theoretical cooling curves with observations. 

Superfluidity of neutrons and/or protons in neutron- 
star cores affects the heat capacity of nucleons and 
reduces neutrino reactions (Urea and nucleon-nucleon 
bremsstrahlung processes) involving superfluid nucleons 
(as reviewed, e.g., by Yakovlev et al. I1999|l . Moreover, su- 
perfluidity initiates an additional neutrino emission mech- 
anism associated with Cooper pairing of nucleons (Flowers 
et al. Q976). All these effects of superfluidity are incorpo- 
rated into our cooling code. 

In our calculations we adopt one model of strong su- 
perfluidity of protons (with the maximum of T cp (p) about 
T™ ax «7x 10 9 K) and several models of moderate super- 
fluidity of neutrons (with T™ ax ~ 6x I0 8 K) in a neutron- 
star core. These models are phenomenological but con- 
sistent with the results of microscopic theories. A pair of 
models: proton superfluidity pi and neutron superfluidty 
ntl is plotted in the left panel of Fig. ^ 

The strong proton superfluidity is required to slow 
down the cooling of low-mass stars, M <^ 1.1 Mq, whose 
central densities are p c ^ 8 x 10 14 g cm -3 . This scenario 
was suggested by Kaminker et al. (2001). In a low-mass 
star, one has T c (p) <; 3 x 10 9 K everywhere in the core. 
The proton superfluidity occurs at the early cooling stage 
(t ^ 1 yr) and suppresses modified Urea processes of neu- 
trino emission as well as neutrino generation in proton- 
proton and proton-neutron collisions. Neutrino emission 
due to Cooper pairing of protons is switched on too early 
and becomes inefficient in middle-aged neutron stars we 
are interested in. In contrast, the adopted neutron super- 
fluidity is too weak in low-mass stars (the left panel of Fig. 
to appear at the neutrino cooling stage. This superflu- 




Fig. 1. Left: Density dependence of critical temperature of model pi for proton superfluidity and model ntl for neutron 
superfluidity in a neutron-star core; vertical dot-and-dash line indicates the central density of a maximum-mass neutron 
star. Right: Observations (Table 1) compared with theoretical cooling curves of eight neutron stars (1-8) with different 
masses. All solid curves refer to neutron stars with model superfluidities from the left panel. The dotted curve 7 is for 
a non-superfluid star. Insert table gives masses and central densities of stars 1-8. 



idity does not suppress the neutrino emission in neutron- 
neutron collisions which becomes the leading mechanism 
of neutrino cooling. It is much weaker than the modified 
Urea process (which would be leading in non-superfluid 
stars). As a consequence, the cooling curves of low-mass 
stars go noticeably higher than the universal cooling curve 
of non-superfluid stars. Actually, these cooling curves also 
merge into one almost universal curve, which is indepen- 
dent of the equation of state in a stellar core and the exact 
behaviour of the T cp (p) (Kaminker et al. 2002). This up- 
per curve 1 allows one to explain observations of neutron 
stars hottest for their age (RX J0822-4300, IE 1207.4- 
5209, PSR B1055-52, RX J0720.4-3125) as cooling low- 
mass neutron stars. 

Now we come to observations of neutron stars cold- 
est for their age (first of all, PSR J0205+6449, the Vela 
pulsar, and Geminga). It has been widely proposed to in- 
terpret these objects as rather massive neutron stars with 
the neutrino emission enhanced by direct Urea process in 
nucleon cores (or by similar processes in pion-condensed, 
kaon-condensed or quark cores) . We will show that coldest 
objects can be explained without invoking these mecha- 
nisms by tuning the model of moderate neutron superflu- 
idity at p ^ 8 x 10 14 g cm" 3 . Let us consider the most 
massive neutron star (1.994 M©, curve 8) in Fig. ^ Its 
central density is higher than the density at which neu- 
tron superfluidity ntl dies out. When the internal tem- 
perature of the star becomes lower than the maximum 
critical temperature of the neutron superfluidity, the neu- 
trino emission due to Cooper pairing of neutrons switches 
on and becomes a powerful neutrino emission mechanism, 



which can be about two orders of magnitude more effi- 
cient than the modified Urea process in a non-superfluid 
star (see Sect.^J). This emission produces enhanced cooling 
(attributed to direct Urea or similar processes in previous 
calculations). The enhancement is not too strong (e.g., the 
direct Urea process in a nucleon stellar core would further 
enhance the neutrino luminosity by about 4-5 orders of 
magnitude). However, even this not too strong enhance- 
ment is sufficient to explain observations of the coldest 
neutron stars (particularly, PSR J0205+6449, the Vela 
and Geminga pulsars). Evidently, all neutron stars with 
p c > 2 x 10 15 g cm~ 3 (in our model) will cool nearly as 
fast as the 1.994 M© star in Fig.Q 

Therefore, we come to three distinct classes of cool- 
ing neutron stars (similar to those described by Kaminker 
et al. 120021 for the case of enhanced cooling due to di- 
rect Urea process). The first class contains low- mass, very 
slowly cooling stars (curve 1 in the right panel of Fig. 
^ . Another class contains high-mass stars with enhanced 
cooling (curve 8). Finally, there is a class of medium-mass 
neutron stars (curves 2-6) which show intermediate cool- 
ing. Their cooling curves fill in the space between the 
upper curve for low-mass stars and the lower curve for 
high-mass stars. These curves explain observations of PSR 
B1706-44, PSR J0538+2817, and RX J1856.4-3754. 

4. Cooper-pairing neutrino emission as a 
fast-cooling agent 

Let us give a simple explanation of the computer results on 
enhanced neutrino emission due to Cooper pairing of neu- 
trons. We start from the expression for the neutrino emis- 
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sivity Qcp due this process (e.g. 
al. I2(JU1|I . It can be written as 

Qcp( P ,T) = q(p,T)F(T), 

where 

q(p,T) « 1.17 x 10 



Eq. (236) in Yakovlev et 40 



(1) 



( m *N 



xTc \M V o/v erg cm s 



PF 

m N c 
l 



(2) 

T = Tg x 10 9 K is the internal stellar temperature, mjy 
is the bare nucleon (N = n or p) mass, m* N is the nu- 
cleon effective mass in dense matter, pp is the nucleon 
Fermi momentum, a at is a dimensionless constant com- 
bined of squared weak-interaction constants of vector and 
axial-vector nucleon currents, M v = 3 is the number of 
neutrino flavors, and F(t) is a function of r = T/T c . The 
constant a at depends on nucleon species and pairing type, 
while F(t) depends on pairing type. We have a n = 4.17 
for the triplet-state neutron pairing under discussion. This 
value can be renormalized by many-body effects (for in- 
stance, the renormalization of the axial-vector constant 
was considered by Carter & Prakash 2002). However, the- 
oretical cooling curves are not too sensitive to the exact 
value of a n , and we use the non-renormalized value. The 
analytic fit expression for F(t) is presented, for instance, 
by Yakovlev et al. |2QQ1J[1. Let us remind that F(t) w 
4.71 (1 — r) just after superfluidity onset (immediately af- 
ter T falls below T c ) and F(t) « 1.27 t~ 6 exp(-2.376/r) 
at t -C 1. Thus, the emissivity Qcp(p,T) is exponentially 
suppressed at T«T C . 

For our qualitative analysis in this section we employ 
the simplest dependence of the neutron critical tempera- 
ture on distance r from the stellar center: 
(r-r m ) 2 



T cn (r) = T a 



1 



(3) 



T™ ax at 



(Ar m ) 2 

at \r — r m \ < Ar m (with the maximum T cl 
r = r ra ), and T cn = at \r - r m \ > Ar m . 

Neglecting, for simplicity, general relativistic effects 
and assuming an isothermal stellar core at a temperature 
T < Tom, the neutrino luminosity Lcp due to Cooper 
pairing of neutrons can be written as 



^CP 



47T 



' Qcp dr. 



(4) 



Here, T\ and r 2 restrict the superfluid layer, where T < 
T cn and the neutrino process in question is allowed. To 
be specific, let us assume that the widest superfluid layer 
(which is realized at T = and extends from r m — Ar m 
to r m + Ar m ) entirely falls in the neutron star core. 

The factor F(t) in the emissivity Qcp, Eq. QJ, is a 
more rapidly varying function of r than q(p, T). Thus we 
can set r = r m and q(p, T) = q(p m , T) (with p m = p(r m )) 
in all functions under the integral but in F(t). A simple 
replacement of integration variable leads then to 

L C p = 8-Krl l Ar m q(p m ,T)T m £(T m ), 



£(t) 



At' F{t') 

r ,3/2/p-Z 



(5) 
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Fig. 2. A sketch of neutrino luminosities produced by 
the modified Urea process (iMurca) and Cooper pairing 
process (Lcp) as well as of the photon luminosity L 7 of 
a neutron star versus internal temperature T for three 
models of neutron superfluidity in the stellar core with 
T™ ax = 10 8 , 3 x 10 8 and 10 9 K. 



where r m = T/T cm . The integration can be done numeri- 
cally; the appropriate analytic fit (for triplet-state neutron 
pairing) is 



a 



x3/2 



3.844(1 
13.99t(1 - r) + 



t) +3.142 t 1 
25 At 2 - 5 (1 - 



((r - 0.2493) 2 + 0.03694) 



0.7 



(7) 



Evidently, the luminosity L^p vanishes in a hot star, 
where T > T cm and neutron superfluidity is absent. 
It switches on as T falls below T cm ; it grows almost 
linearly while T decreases to ~ 0.8 T cm ; afterwards, it 
reaches maximum at T = 0.792 T cm (with t£(t) = 
0.792^(0.792) = 0.481) and then decreases. At the increas- 
ing and maximum-luminosity stage, Lop is collected from 
a superfluid spherical stellar layer in the vicinity of the 
maximum critical temperature, r w r m . This creates a 
splash of neutrino emission associated with Cooper pair- 
ing of neutrons. 

For typical values of the parameters, the maximum 
value of Lcp can be one-two orders of magnitude higher 
than the neutrino luminosity LMurca of a non-superfluid 
star (with forbidden direct Urea process). This is demon- 
strated in Fig. |2 using a toy model of cooling neutron 
stars described by Yakovlev & Haensel i2003J) - there is 
no need to employ accurate models in this section. The 
parameters of the neutron-star model presented at the 
figure are: M = 1.16 Mq, R = 12 km, p c = 8 x 10 14 
g cm -3 , r m = Ar m = 5 km. The three superfluidity 
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models (T cn (r)) are selfsimilar and differ by the values 
of T cm = 10 8 , 3 x 10 8 , and 10 9 K. Three solid lines ex- 
hibit the Cooper-pairing neutrino luminosity Lcp calcu- 
lated from Eqs. Q-Q for three models of neutron super- 
fluidity. Since L M mca ^ 8 an d L™p x 00 ■, the Cooper- 
pairing luminosity is more competitive at weaker super- 
fluidity (lower Tcm). However, at T cm < 2 x 10 s K this 
luminosity becomes lower than the photon thermal lumi- 
nosity of the star (Fig. [21 which makes it insignificant 
for stellar cooling. It is worth to notice that, for realistic 
parameters, Lcp is much smaller than the neutrino lumi- 
nosity due to the direct Urea process in a non-superfluid 
star (if the direct Urea process is open). 

The decreasing part of Lcp (T) is even more fascinat- 
ing. We have i{r) ~ 3.84 as r — •> 0, resulting in the scaling 
relation 

Lcp oc Ar m T s /T cm , (8) 

which becomes sufficiently accurate at T < 0.6 T cm . This 
neutrino emission is actually produced from two thin 
spherical shells (near r = r% and r = r%), where T is 
just below T c (r). The widths of these shells are propor- 
tional to T, which explains the power-law T 8 (instead 
of the exponential decrease of the emissivity Qcp(p,T) 
in a local element of superfluid matter). Therefore, the 
decreasing part of the Cooper-pairing neutrino luminos- 
ity has the same temperature dependence as all slow 
neutrino emission mechanisms (modified Urea, nucleon- 
nucleon bremsstrahlung) in non-superfluid cores of stars 
with forbidden direct Urea processes. In other words, the 
superfluidity suppresses the neutrino emission available in 
non-superfluid stars but initiates the Cooper-pairing neu- 
trino emission in such a way that it acts as a new nonsup- 
pressed neutrino cooling mechanism. Moreover, the new 
emission can be more intense than that in a non-superfluid 
star and provide enhanced cooling. This important feature 
appears in realistic models of cooling neutron stars with 
density dependent critical temperatures T c (p) (and does 
not appear in the models with density- independent T c ). 
In particular, it implies that once Lcp takes on leadership 
in competition with LMurca just after the superfluidity on- 
set, it will not lose it during subsequent evolution (espe- 
cially because LMurca is actually suppressed by superflu- 
idity, which is not taken into account in Fig. (2J. This is 
clearly seen from Fig. 

Let us add that at T < T cm we can obtain a bet- 
ter formula for Lcp than Eq. J3J, without employing the 
specific T c (r) profile, Eq. (JSJ. It is sufficient to start from 
Eq. Q and notice that the main contribution into Lcp 
comes from two thin shells, at r « r\ and r ~ r%, where 
T c (r) « T. In each shell, the gradient D = dT c (r)/dr can 
be taken constant. Then we get 

L CP = 8^ [r? #i q( P i,T) + r\ H 2 q(p 2 ,T)] £(0), (9) 

where H\ — T/\D\\ and H2 = T/\D2\ are characteris- 
tic widths of our shells, p\ = p{r\), pi = p{r<z), and 
^(0) = 3.84. Strictly speaking, n, r-x, pi, p%, D±, and 



L> 2 depend slightly on T, but this dependence can be re- 
garded as parametric. It is easy to verify that if T c (r) is 
given by Eq. J2J at T <C T cm and Ar <C ?" m , our new ex- 
pression for Lcp coincides with Eq. JSJ. Equation JSJ is 
expected to be useful just after the superfluidity onset, at 
0.6 T cm 5s T < T cm (where the parabolic T c (r) dependence 
may be a good approximation), while Eq. is more ex- 
act at lower T. Both equations enable one to incorporate 
the Cooper-pairing neutrino emission in simplified cooling 
models (like a toy model of Yakovlev & Haenscl 2003 ) , use- 
ful for understanding main features of neutron star cooling 
without complicated cooling codes. 

The above analysis is valid as long as T c (p) vanishes in 
the stellar interior. If not, there is a minimum value T c mm 
of T c {p), and Lcp will become exponentially suppressed 
at T < T c min . 

5. Testing the cooling scenario and discussion 

After clarifying the efficiency of the Cooper-pairing neu- 
trino emission let us return to the cooling scenario de- 
scribed in Sect. 01 As we have already mentioned, the sce- 
nario is rather insensitive to a specific model of proton 
superfluidity (required to raise the surface temperature of 
low-mass stars for explaining observations of the sources 
hottest for their ages). The only serious constraint on the 
proton pairing is that T cp (p) should be high (<; 3 x 10 9 
K) in the cores of low-mass stars. 

However, the constraints on the neutron critical tem- 
perature T cn (p) in a stellar core should be really strong. 
This is illustrated in Fig. [31 The left panel displays the crit- 
ical temperatures of our proton superfluidity model (pi) 
and five neutron superfluidity models (ntl-nt5), includ- 
ing our basic model ntl used in Sect. 03 The right panel 
shows cooling curves of a low-mass (1.111 Mq) star and 
a high- mass (1.994 M ) star. Any curve is calculated for 
model pi of the proton superfluidity and one model of the 
neutron superfluidity from the left panel of Fig. Any 
observational point between an upper curve and a lower 
curve can be explained by a given superfluid model. The 
constraints on the neutron superfluidity are as follows. 

(1) The neutron superfluidity should be weak in low- 
mass stars. In our case (for the equation of state of 
Douchin & HaenselEEHJl this means that T cn (p) 5s 2 x 10 8 
K at p 5s 8x 10 14 g cm" 3 . Under this condition the neutron 
superfluidity does not affect the cooling (at least at the 
neutrino cooling stage) of low-mass stars [M 5s 1.1 M©) 
and does not violate our interpretation of the sources 
hottest for their age (first of all, RX J0822-4300 and PSR 
B1055-52). Accordingly, all five cooling curves (for su- 
perfluids ntl— nt5) of low-mass stars merge in one upper 
(solid) cooling curve in Fig. |3 The only exclusion is pro- 
vided by model nt2 with highest pre-peak T cn (p) among 
models ntl-nt5. In a low-mass star this superfluidity oc- 
curs at t ^ 300 kyr. The Cooper-pairing neutrino emission 
and reduced heat capacity of neutrons noticeably acceler- 
ate the cooling at this late stage (the upper short-dashed 
curve). 




Fig. 3. Left: One model pi for proton superfluidity and five models ntl-nt5 for neutron superfluidity in a neutron-star 
core. Right: Cooling curves of low-mass (1.111 Mq) and high-mass (1.994 Mq) stars with model pi proton superfluidity 
and one of the models of neutron superfluidity from the left panel, compared with observations. Cooling of the low- 
mass star is insensitive to selected models of neutron superfluidity (except for model nt2 at t > 300 kyr) . Insert shows 
the comparison of cooling curves of the high-mass star with observations of the Vela pulsar in more details. 



(2) The neutron superfluidity should be moderate at 
p ^ 10 15 g cm -3 , with the peak maximum T™ ax ~6x 10 s 
K (model ntl in Fig. the solid curve). In this case it 
switches on just in time to initiate the enhanced cooling 
in a high- mass star. Its level is sufficient to explain obser- 
vations of neutron stars coldest for their ages (first of all, 
PSR J0205+6449 and the Vela pulsar). The asymptotic 
neutrino-cooling regime given by the scaling expression (JSJ 
is realized at t jsj (1 — 10) kyr. If T™ ax were slightly higher 
than 6 x 10 s K (model nt3, T™ ax = 8 x 10 8 K, the dotted 
curve), the Cooper-pairing neutrino emission will start op- 
erating in a younger massive star but becomes less efficient 
at t ~ 10 kyr, which is less favorable for explaining the 
observations of the Vela pulsar. This cooling behaviour 
is naturally explained by the scaling (JjJJ. If T™ ax were 
slightly lower than 6 x 10 8 K (model nt4, T™ ax = 4 x 10 8 
K, the long-dashed curve), the Cooper-pairing neutrino 
emission will start operating too late which would violate 
the interpretation of the observations of PSR J0205+6449. 

(3) The results are also sensitive to the width of the 
peak of the T cn (p) curve. For instance, retaining the peak 
maximum of 6 x 10 8 K but making the peak narrower 
(model nt5, the dot-dashed curve) will reduce the neutrino 
emissivity due to neutron pairing, raise the temperature 
of the massive star and complicate the interpretation of 
the Vela pulsar (again, in agreement with the scaling (JSJ). 
However, the cooling curves are rather insensitive to the 
exact position of the T cn (p) maximum. We can slightly 
shift the maximum to higher or lower p (confining the 
peak within the kernel of a massive star) but these shifts 
will not change the cooling curves of massive stars (such 



tests are not shown in Fig. |3J). However, the shift of the 
maximum to p ^ 8 x 10 14 g cm -3 would cause the en- 
hanced cooling of low-mass stars. The cooling curves of 
low-mass stars would become close to those of high-mass 
stars which would violate the interpretation of neutron 
stars hottest for their ages (see item (i)). 

This discussion shows that the cooling curve of a mas- 
sive neutron star implying model ntl of neutron superflu- 
idity is close to the lowest cooling curve (in the scenario, 
where the cooling is enhanced by Cooper-pairing neu- 
trino emission). Observations of cold neutron stars, PSR 
J0205+6449 and the Vela pulsar, provide excellent tests 
for this scenario. If these pulsars were noticeably colder 
we would be unable to explain them within our scheme. 
Notice, that the upper limit of the surface temperature of 
PSR J0205+6449 was inferred from observations (Slane et 
al. I2002fl using the blackbody spectrum of surface emis- 
sion. If this pulsar has a hydrogen atmosphere, the upper 
limit on T s °° could be expected to be about twice lower 
than for the blackbody case. In that case we would be un- 
able to explain this source within the proposed scenario. 

Although we have used one equation of state of dense 
matter (Douchin & Haensel 2 001(1 we would obtain simi- 
lar results for other equations of state which forbid direct 
Urea processes (and other similar processes of fast neu- 
trino cooling) in neutron star cores. Taking different equa- 
tions of state would lead to attributing different masses to 
the same sources (Fig. ^ ; similar problem has been dis- 
cussed by Kaminker et al. (2002). 

In addition, we could take an equation of state in the 
stellar core which opens direct Urea process at highest 
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densities (in the central kernels of most massive stable 
neutron stars; similar to the equation of state of Akmal 
& Pandharipande lT997[) . Applying the same model of nu- 
cleoli superfluidty as in Fig. ^ we would get five types 
of cooling neutron stars (instead of three). Three types 
would be the same as those mentioned in Sect.|31 low-mass, 
very slowly cooling stars; massive stars whose cooling is 
enhanced by Cooper-pairing neutrino emission; and of 
medium-mass stars whose cooling is intermediate. In ad- 
dition, we would have: most massive neutron stars demon- 
strating a very fast cooling via the direct Urea process; and 
stars whose cooling is intermediate between that enhanced 
by the Cooper-pairing neutrino emission and by the direct 
Urea process. The transition from the Cooper-pairing neu- 
trino cooling to direct-Urea cooling with increasing mass 
M will be very sharp and the number of intermediate- 
cooling sources will be small. The maximum-mass neutron 
stars would be extremely cold (T s °° ~ 2 x 10 5 K at t ~ 10 
kyr), about the same as discussed, e.g., by Kaminker et al. 
( 2002J. A discovery of such stars would definitely indicate 
the operation of the direct Urea process in their cores. An 
indirect evidence of their existence is provided by a non- 
detection of neutron stars in some supernova remnants 
(Kaplan et al-EH®. 

Note that the cooling of neutron stars can also be af- 
fected by the singlet-state superfluidity of neutrons in in- 
ner stellar crusts, by the presence of surface layers of light 
(accreted) elements, and by stellar magnetic fields (e.g., 
Potekhin et al. EHH Geppert et al. 1271714) . These effects 
can be especially important in low-mass stars. We have 
neglected them in the present paper since we have mainly 
focused on enhanced cooling of massive stars but we will 
consider them in a future publication. 

6. Conclusions 

We have proposed a new scenario of cooling of isolated 
neutron stars. We have shown that the present observa- 
tional data on thermal emission from isolated middle-aged 
neutron stars can be explained assuming that neutron 
star cores are composed of neutrons, protons and elec- 
trons (and possibly muons) with forbidden direct Urea 
process of neutrino emission. In our scenario, an enhanced 
neutrino emission, which is required for interpretation of 
neutron stars coldest for their age, is provided by neu- 
trino process associated with Cooper pairing of neutrons. 
We have shown that the neutrino luminosity due to this 
process (at internal temperatures T ^ 0.6T™ ax ) behaves 
as T 8 . In this way it "mimics" the neutrino luminosity 
produced either by modified Urea processes or nucleon- 
nucleon bremsstrahlung processes in non-superfluid stars, 
but it can be one- two orders of magnitude higher. The pro- 
posed cooling scenario imposes very stringent constraints 
on the density dependence of neutron-pairing temperature 
T cn (p). The constraints result from the comparison of the 
cooling theory with two most important "testing sources" , 
PSR J0205+6449 and the Vela pulsar (Sect.®. This sce- 
nario is the first one in which a moderate superfluidity 



and associated neutrino emission are helpful for explain- 
ing the data (cf. with previous cooling scenarios, where a 
moderate superfluidity has violated interpretation of ob- 
servations, e.g., Kaminker et al. [2002). 

Our interpretation implies the presence of a strong pro- 
ton superfluidity and a moderate neutron superfluidity in 
neutron star cores (Sect.®. We need the proton superflu- 
idity to explain observations of neutron stars hottest for 
their age, and the neutron superfluidity to explain obser- 
vations of stars coldest for their age. However, as has been 
demonstrated by Gusakov et al. (2004), cooling curves 
are not too sensitive to exchanging neutron and proton 
superfluidities (T cp (p) # T cn (p)) in neutron-star cores. 
Therefore, we would also be able to explain observational 
data in the scenario with a strong neutron superfluidity 
and a moderate proton superfluidity in stellar cores. 

We need a strong superfluidity to suppress modified 
Urea process in low-mass stars, rise the surface temper- 
ature of these stars and explain observations of neutron 
stars hotter for their age. In fact, we can rise the temper- 
ature of low-mass middle-aged neutron stars by assuming 
the presence of surface layers of light (accreted) elements. 
The mass of light elements may decrease with time, e.g., 
due to diffusive nuclear burning (Chang & Bildsten[5003), 
which opens additional freedom to regulate the cooling. In 
this way, the presence of a strong (proton or neutron) su- 
perfluidity in a neutron star core is not vitally important 
for our interpretation. We will show this in a future publi- 
cation. However, the presence of a moderate superfluidity 
(of neutrons or protons) with a tuned density dependence 
of critical temperature (Sect. ® is crucial for this sce- 
nario, where this tuned dependence is combined with the 
remarkable simplicity of the equation of state of neutron- 
star cores (nucleon composition with forbidden direct Urea 
process). We hope that this scenario can be taken into 
consideration along with many other scenarios (reviewed 
or proposed, e.g., by Page I199"5al 11998b! Tsuruta et al. 
12717121 Khodel et al. 1271711 Blaschke et al. 1271711 Yakovlev & 
Pethick 2004j) . The correct scenario should be selected in 
future observations of neutron stars combined with new 
advanced theoretical results. 

After this paper was prepared for submission we be- 
came aware of the paper of Page et al. ( 200||). These au- 
thors give a detailed consideration of enhanced cooling 
via neutrino emission due to Cooper pairing of neutrons 
in neutron-star cores composed of nucleons with forbid- 
den direct Urea process. The idea to enhance the cooling 
by Cooper-pairing neutrino emission is the same as in our 
paper, but its realization is different. Particularly, Page 
et al. H2()()4|l use a set of superfluidity models obtained 
from microscopic theories. Their main models for neutron 
superfluidity in a stellar core (for, instance, model (a) in 
their Fig. 9) have too high peak temperatures T c ™ ax ^ 10 9 
K and too high T cn (p) at the pre-peak densities to explain 
the observations of PSR J0205+6449 and the Vela pulsar 
and to obtain a pronounced dependence of cooling curves 
on neutron star mass. In contrast, our T cn (p) models are 
phenomenological but tuning them we obtain a noticeable 
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dependence of the cooling on M. It enables us to explain 
all the data by one model of nucleon superfluidity (even 
neglecting the effect of accreted envelopes). 
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